arXiv:1501.01092v2 [math.AP] 7 May 2015 


DUALITY THEORY AND OPTIMAL TRANSPORT FOR SAND 
PILES GROWING IN A SILOS 


LUIGI DE PASCALE AND CHLOE JIMENEZ 

Abstract. We prove existence and uniqueness of solutions for a system of 
PDEs which describes the growth of a sandpile in a silos with flat bottom 
under the action of a vertical, measure source. The tools we use are a discrete 
approximation of the source and the duality theory for optimal transport (or 
Monge-Kantorovich) problems. 


1. Introduction 


Let be a bounded, convex open subset of with 2 < d G N, let : dQ —>■ M"*" 
be a bounded lower semi-continuous function, and let / G L°°(0, T, AI(r2)) with 
0 < /• We denote by AI(r2) the set of Borel measures on hi with finite total 
variation. When / is time constant we will prove existence of solutions {u, /i, u) of 
problem 


(PDE) dtu — div(D^jUpf) = / — z/ in R'^xjO, T[, 

(CV Constraints I ^ nx]0,T[, 

(B): Boundary conditions (<((1) 0 <“U‘) < 9(d in afix]0 r| 

^ ^ 1(52) u{x,t) = g{x) ut-&.e.{x,t) e D.x]0,T[, 


(I); Initial conditions 'u(-,0) = 0 in hi, 
with the regularity fll.2|]ll : 


( 1 . 1 ) 


u G L°°(0,T,IU^’°°(fi)), dtU G L°°(0,T,M(fl)), 
/i G L°°(0,r, Al+(fi)), z/ G L°°{0,T,M+{dQ)). 
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{x,t) for /xt-a.e. x, a.e. t. 
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A function u as above is, in particular, in C([0, T]; which gives sense to 

the initial data f fll.lD .fllb 

In these equations Du denotes the spatial part of the derivative while dfU is 
the time derivative of u. The divergence is only spatial and is intended in the 
distributional sense. The meaning of will be given later and should be seen 
as the part of Du which is relevant for the measure fi. 

As usually the equality given by (PDE) must be intended as 

^ [ u{-,t)ip{-) dx+ [ D^,^u- Dip diit= [ (pdft- [ (f dut, mV'{0,T) 
jR'* J Jn Jan 

for all if G D{W^). Here D{A) denotes the space of smooth functions compactly 
supported inside A. 

The function u will be proved to be unique, we will show that the problem above 
is equivalent to a variational inequality in the spirit of the original paper j^H]- One 
of the main tools will be the duality theory for optimal transport problems. This 
approach will allow us to study some form of uniqueness and mild regularity for 
p. More precisely, we will prove the two following results: 


Theorem 1.1. Let f G L°°(0, T, Al(r2)) and {u,p,iy) satisfying U.^) with the 
boundary condition U.li fBD). and assume u{-,t) ^ Cipi{Q) a.e. t. 

Then [u, p, u) is a solution of U-l\) if and only if (I) is satisfied and ft — dtu{-,t) ^ 
dIoo{u{-,t))a-e. t g]0,T[. 

This last condition means thatu is a solution, a.e. t, of the following maximization 
problem: 


ma.x{{f — dtu{-,t),v) : v E Cipi{D), 0 < v{x) < g{x) on dfl}. (1.3) 

Moreover M.Al is the dual formulation of the mass transport problem and to 
any optimal choice of u corresponds a unigue optimal p defined by: 
r /•! 


{Tt,T) 


V3((l - s)x + sy)\y - x\ ds d'yt{x,y) Mp G Cfe(n) 


Jn Jo 

where (together with Ut) is any solution of Ii4-4\ l- 

Theorem 1.2. Let f G L°°{0,T, Ai{D)) be constant in time. Then M.l\) admits 
a solution {u, p, v) satisfying U.Sf) . Moreover u is unigue and dtu G L^(r2x]0, T[). 


The results above will be proved by approximating / by a finite number of point 
sources. 

The differential system (HHi has been proposed several years ago in m to 
describe the growth of a sandpile on a bounded table under the action of a vertical 
source here modeled by /. In the model the sandpile is described as composed 
by an underlying standing layer, here modeled by u and a rolling layer which is 
here modeled by p. The material rolls downhill only when the standing layer 

^see, for example, [55] page 54 
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reaches a critical slope which is characteristic of the material. This critical slope, 
here, is normalised to 1 and the conditions (C) gives account of this aspect of the 
behaviour. Together with the boundary conditions {BC) the system models the 
growth of the pile inside a silos with wall on dfl of height g. At some point the 
sandpile will reach the top of the wall and the sand will start to fall out. The 
measure v describes where this will happen and how much sand will fall out from 
each point. So at the beginning we expect Pt to be 0 while after some time the 
sandpile will stabilise and dtu will become 0 and v will have the same mass of /. 

Allowing V to have sign (which is equivalent to assume that an additional source 
of sand may appear on the boundary) could bring to a loss of uniqueness as shown 
in the following example. 


Example 1.3. Take d = 1, fl =]0,1[, / = 0, T = 1 and g constantly equal to 1. 
Obviously (u, p, u) = (0, 0, 0) is a solution, but we may also take for instance: 



u{x,t) = {t - x)l{^^^t):x<t}{x,t), Pt = 


t rl 


(p(x, t) dgt{x) = 



(p((l — s)x)x dsdx 


0 ^0 


—tho, 

V^ea(]0,l[). 


The choice of / in the space of measures aims to model situations in which 
the source has dimensions much smaller than the sensitivity of the measure in¬ 
struments. A possible example from daily life is a hourglass where sometimes the 
passage for the sand is small at the limit of the imperceptible. 

A different approach to this problem is currently pursued by other authors |18) . 


Remark 1.4. When / is a measure it may happens that g is not better than a 
measure as shown in several examples in 

As every model, this one is well suited for some situations and it fails in others. 
An accessible description of several models (included the one we consider), may 
be found in [27] together with several more references. 

Most of the literature is concerned with the Dirichlet case (also known as table 
problem) with a source / G U{VL). In that case the system describes a sandpile 
growing on a table without walls. At some moment the pile reaches the boundary 
of the table and the sand start to fall out stabilizing the pile. The standard 
approach consisted in proving that in a regular setting the system is equivalent to 
a variational inequality which may be written as 

f - dtU e dIoo{u) (1.4) 

and then proving the existence and uniqueness of a solution u of fll.4p for a wide 
class of /. For example, this has been done by Prigozhin in [29] where existence 
and uniqueness of u (for the table problem without walls) is proved under the 
assumption / G (L^(0, T, kF^’^(n))' which is wider than the space we consider. 
Nevertheless, the meaning of g is less explicit in that work. Making the link with 
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optimal transport, and making nse of the dnality theory, allows ns to give sense 
to /i when / is a measnre. 

A similar model for a sandpile growing nnder the action of a finite nnmber 
of sonrces was proposed earlier in [3] and in several nnpnblished notes by the 
same anthor. The ODE arising from that approach was stndied in [1] and in the 
same paper approximations as p —>■ oo of the problem were also introdnced. 
The same kind of approximations are nsed in [23] to stndy the problem for more 
general sonrces and to establish a first relationship with the optimal transportation 
problem. 

For the table problem a setting similar to the one of this paper is nsed in 
PH1I25] for a sonrce / G T^(0, T,A4{Q)). The same papers contain some theory for 
nnmerical approximations as well as nnmerical simnlations. One of the difhcnlties 
one has to face when / is a measnre (in the space variable) is dne to the lack 
of regnlarity of p. The approach of N. Igbida avoids this problem by working 
with the flnx and introdncing some weak formnlation of fll.ip . To deal with this 
lack of regnlarity of p, we nse the tangential calcnlns with respect to a measnre 
introdnced by Bonchitte, Bnttazzo and Seppecher in j?]. 

Convergence toward eqnilibrinm in the table problem is stndied in [T2] • Then, 
stationary solntions are stndied, with or withont explicit mention of the sand piles 
model, in [SI El m EDI Ell EZl ESI ES]- There is still mnch work to do toward a 
complete nnderstanding of the time of convergence to eqnilibrinm when the sonrce 
is not too regnlar or controlled from below. The literatnre on the silos problem is 
not so rich. Some partial resnlt together with nnmerics is contained in |13] . 

Here we choose to start from an empty silos. The theory wonld be the same (np 
to technical details) if we assnme any initial condition uq which satisfies \Duo\ < 
1. If \Duq\ = 1 on a set of positive measnre then one enters the domain of 
collapsing sand piles which is different from the one we are considering and is rich 
of interesting problems (we snggest to start from [23] and to follow with the papers 
in which that paper is cited). 

Finally a remark on the convexity of D. It is clear that the problem wonld be 
interesting also in non-convex domains. Some resnlts for stationary (eqnilibrinm) 
solntions is contained in the recent paper [16] which also consider some anisotropic 
generalization. Here the assnmption of convexity is crncial in section 5 where we 
use shortest line connecting internal points of D to some point on dQ determined 
by the values of g. 

Notations : The Euclidian norm on will be denoted by | ■ |. For any A subset 
of ]0, T[ or T[, we denote by: 

• Ai{A) the bounded measures supported in A, M.^[A) being the subset of 
Ai(A) of non-negative, bounded measures, 

• For any vectorial measure a G Ad(A,M"), we will denote by |cr| the total 
variation measure associated to a, 
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• C{A) the set of continuous functions on A while Cb{A) denotes the set 
of continuous bounded functions on A equipped with the inhnity norm 
II • I loo, the smooth compactly supported functions on A will be denoted by 
V{A) ■.= C^{A). 

For any functional space T, we will use T' for the topological dual of T. Slightly 
abusing notations, we will write (p, for (p, (Ct{A)y, cua) for any (p, (p) e 

M{A)xCb{A). 

2. Tangent space to a measure and integration by parts 


When u is not regular enough the product of Du and a measure p does not 
make, a priori, any sense. Indeed Du may be not dehned on a set which has 
positive p measure. Here we report few useful notations and results from |7j. This 
will give sense to D^^u appearing previously. 

Let T] G We can set as in [7]: 

:= {'ll; e : div(i/;r7) G 

Tri{x) ■= T] — ess U {'0(a;) : -0 G Xn}, 

where the divergence is intended in the sense of measures. More precisely "0 G 
is in X^^ iff there exists a constant K such that: 


Dip{x) ■ 'ijj{x)d'q{x) < iL||(p||oo V(p G ViW^ 


For ip G the tangential gradient to ?7 at x of (p is dehned as: 

D^iplx) := Prj{x){Dip{x)) //-a.e.x, with Prf{x) := orthogonal projector on T^{x). 

As shown in [7], the operator u G i—)■ D^u G can be extended by 

setting: 

w =: Dr/V -v^ {3vn G — )■ n uniformly, Dr/Vn ^ tc in (T^)'^). 

The tangential Sobolev space is then dehne as the domain of D^. By dehnition 
any vector held in belongs to the dual H~^ and then the following integration 
by parts formula holds 


DrfU{x) ■'ijj^x) d'q{x) = for all u & and i/’ G X^. (2.1) 


Example 2.1. Assume 77 G and let u G Denote by v any contin¬ 

uous, compactly supported extension of u to It is easily seen that v belongs to 
and that DnV will be the same for any other compactly supported extension of 
u. We denote by D^^u := Dj^v. We have: \Dr^u{x)\ < ||Dm||oo 77 -a.e.x. Moreover, 
m rewrites as: 



D riU{x) ■ '0(x) d'q{x) 


{-d\v{'il;vi),u)M{n),c^{Q.) for all G X, 
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Remark 2.2. In order to build some notion of space-time tangential gradient, we 
could have chosen some notion of space-time tangential gradient, to set for any 
p e L°°(]0,r[,>l(M'^)); 

X, = {^^eLl{W^x]0,T[Y : div(i/./i)eL°°(0,T,7U(R''))}. 

Nevertheless, the result will be exactly the same. Indeed, for any 'ip G L^(M'^x]0, T[)‘^, 
the following equivalence holds; 

p) ^ for a.e. t g]0, T[. 

Let us prove this equivalence. Assume that t) G for a.e.f g]0, T[ and take 
(f G P(R'^), h G 'D(]0,T[) then; 


D{hLp){x,t) ■ = / h{t) / Dip{x) ■'ilj{x,t)d^t{x)dt 


<K h{t)\\ip{-)\\^dt<K \\h{t) Lp{-)\\^ dt. 


This means exactly that pj G XtJ.- The other implication is straightforward. 


3. Duality and optimal transport 

The results contained in j5] (see also ra) suggests to consider the maximization 
problem (13.1|) dehned below. Its link with our problem will appear clearly in 
Theorem 13.21 The set Cipi { fl ) is defined by 

Cipi{Q) := {n G Cip{R‘^) : v{y) — v{x) < |x — i/| in D x D} . 

Following [5], we prove Proposition 13.11 Theorem 13.21 and Proposition 13.31 

Proposition 3.1. Let u G Cipi{Q) and p G A1(D). Then the following extremal 
values coincide: 

max{(p,n);vj(o)A(u) : v e Cip^{Vt), Q < v < g on dVt) (3.1) 

inf < / d \ a \ + / gdp '^ ; —diver = p — y inRf'X . (3.2) 

o-GAJ(0,R'*),i^GAJ(aO) iJq Jq^ J 

Note that any sequence (n„)„ of admissible applications for (13.ip is uniformly 
bounded on cID and as it is in £zp^(f2), it is uniformly bounded in D. The existence 
of a maximizer is then easily proved. 

Proof. Let us introduce for any {p,q) G C{dfl)‘^ the following perturbation func¬ 
tional; 


H{p,q) :=-sup{{p,v)Min),Cb(^) ■ v e Cipi{Q), v + p < g, 0 < v + q on dQ}. 
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Now compute the Fenchel transform of if on a couple (p*, q*) G M.{dVtY-. 

= sup {{p\p)M{dn),c{dQ) + {q*, q)M{dn),C{dn) - H{p,q)} 

p,q&C{dfl) 

= sup sup {{p*,p) + {q*, q) + {p,v) : v E Cip^{fl), —q<v, p < g — v on 917} . 

p,q V 

= sup < sup{(p*,p) + (g*, (?) : —q<v,p<g — vond^l} + {p,v): v E Cipi{Vt) 

V I p,q 

sup^^piip-q*,v) + {p*,p) : V E Cip^{n), p < g - v} if g* < 0; 

+CX3 elsewhere. 

It is easily seen that H*{p*, q*) ^ +cxo only if p* > 0. Let us assume this condition 
is satisfied. We have: 

H*{p\q*) < {p\ g) + sup {{p - p* -q\v) : v E Cip^{n)} . 

V 

Let us show the opposite inequality. Let {gn)n a sequence in C{dfl) converging to 
g at any point of 917 with < g. Then, by taking p = gn — v we have: 



H*{p*,q*)> lim / gn{x) dp*(x) + sup {{p — p* — q*, v) : n G £lpi(17)} . 


n—>-+00 


fdn 


Finally the equality follows using Fatou’s lemma and we have: 

H*(n* n*) = I {{P-P* -q\v) : v E £lpi(17)} if p* > 0 and q* < 0; 

^ ^ } +CXD elsewhere. 

And, by standard duality (see 0): 

(p*,p) +inf<,g^(f^_] 


H*{p*,q*) = 


1 [ dW\ ■ 

-diver = p - p* - q* in R‘^\ 

Un 

J 


if p* > 0 and q* < 0; 
+00 elsewhere. 


It can be easily proved that H is convex. Let us check that it is l.s.c. Let 
iPnidn) £ C(917)^ converging uniformly to (p, g) G C(917)^. For any e > 0 take 
Vn E CipiiVL) £-optimal for H{pn,qn) that is such that H{pn,qn) > —{p,Vn) — s. 
Possibly extracting a subsequence, we may assume: 

lim inf iL(p„, = lim H (p„, qn ). 

n—>-+oo n^+oo 

From the regularity of Vn and the bounds on 917, we get: 

Vn{x) < sup{|p - xl + g{y) -Pn(p)} < diam(17) + \\g - PnWoo 

yGdQ 


Vn{x) > inf {-\y - x\ - qn{y)} > -diam(17) - ||g„ 

yGofl 


LUIGI DE PASCALE AND CHLOE JIMENEZ 


As {vn)n is an equicontinuous and bonnded seqnence, by Ascoli, a snbseqnence of 
{.^nk)k of (Pn)n conveiges nnifornily to some v G Cip^iVL) admissible for H{p,q). 
Then: 


lirnini H{pn,qn) > hm -{p,Vn^) - e =-{p,v) - e > H{p,q) - e. 

n—>-+oo fc—>-+oo 

By sending e to 0, we get the lower semi-continnity of H. The resnlt then follows 
from the eqnality i^(0,0) = □ 

Theorem 3.2. Let p G Ad(f2). 

(i) Assume that (u,p,iy) G Cipi(i}) x x A4.{dVt) is a solution of 


— div{pD^u) = p — v in 
( n < 

with 


\D^u{x)\ = 1 p-a.e.x, 


Q < u < g on dQ, 

u{x) = 0 i'~-a.e.x, u{x) = g{x) y^-a.e.x. 

Then u G a,rgmax l\d.l\} and, setting a = D^up, we have (a, v) G ar amin {\?).2\\ . 


(3.3) 

(3.4) 


(ii) Moreover, if u E Cipi{fl) and (cr, u) G x AA{dfl) are optimal solutions 

of \3.1\) and \3.2t) then setting p = \a\, we have: 


a = Dfj^up and l-D^nl = 1 p — a.e, 
u = 0 n~-a.e.x, u{x) = g{x) v^-a.e.x. 


(3.5) 


Proof, (i) With the assnmptions above, u and {a = D^up, v) are admissible for 
fl3.ip and fl3.2p . An integrating by parts (see fl2.1|B leads to: 

(p, w)aj(o)A( o) = (-div(/iT>^n) + P, C( 0 ) = j\D^u\^dp+j u du 

Jn Jan 


/ d\a\ + g du^. 
In Jan 


By Proposition 13.11 this implies the resnlt. 
(ii) By Proposition 13.11 


this implies 


and 


(P, w)aj(o) A(o) = dp+ g{x) du^{x) 
Jn Jan 


(-div(a) + = dp + g{x) du+{x), 

Jn Jan 


(-div((7),M)^pc{o) - dp= g{x)du+{x)- / u{x) du{x). 

Jn Jan Jan 
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Finally, integrating again by parts (see again (12.Ih ): 

f {Dfj^u{x) ■—{x)— 1) diJ,{x) = [ {g{x) — u{x))di'~^{x)+ f u{x)du~{x). (3.6) 

Jq dfx Jqq Jqq 

To conclude, we notice that on the one hand as u is in Cipi{Q) 

dcT 

D^u[x) ■ ^(x) — 1 < 0, p-a.e.x 
on the other hand, on dVL 

g{x) — u{x) > 0 and u{x) > 0 . 

This inequalities, combined with fl3.6p give (13.51) . □ 

The following proposition shows that (13. ip can be seen as a variant of the dual 
formulation of the classical Monge mass transportation problem (see [T] or |31) for 
instance): 

Proposition 3.3. Let p G A^(r2). Then the extremal value US. ill coincides with 
the following extremal one: 

_nnn < / \x — y\d'y{x,y) + / gdn'^ \ Trl'y = p~ + , tt ?7 = p"*" + 

■yeM+{nxn),ueMidn) {JnxTi Jan 


(3 


This is a consequence of classic duality in the theory of optimal transport 
(see for instance m) which implies the following equality: 

min (l3.7p = min {Wi{p~ + , p'^ + i'~) + / a di/~^} = inf (l3.2p . (3.8) 

ueM(dn) Jqq 

Where ITi is the Wasserstein distance: 

Wi{p~ + p+ + u~) : = 


mm _ < / lx - yld'y^x, y) : vrK = p + i/+, 71^7 = + 

7eA(+(OxO) IJfixn 


V 


We have the following result: 

Lemma 3.4. (i) Let ( 7 , n) G Wl(r2 x hi) x Ai(dQ) be a solution of LS. 71) . then 

it exists a unique a G Wl(r2,R‘^) such that [a, u) is a solution of 

(o', 0)x(o,Rb — / / (l){{l-s)x + sy)-{y-x) dsdj{x,y) V0 G Cb(f2, M'^). 

702 Jo 

On the contrary, if (cr, v) is a solution of l[3.2\} . it exists 7 such that ( 7 , v) 
is optimal for ^3. 71) and the previous expression of a holds. 
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(ii) Let u optimal for LS. il) and a as in (i). Then u is differentiable \a\-almost 
everywhere and: 

D\a-\u{x) = Du{x) \a\-a.e. x. 

(iii) Let ( 7 , u, u) admissible for LS. 71) and LS. il) . then ( 7 , p, u) are optimal if and 
only if: 

( u{y) - u{x) = \x - y\ 7 - a.e. (x, ?/) 

^ u{x) = 0 v~-a.e.x, u{x) = g{x) u^-a.e.x. 

Proof. Point (i), using fl3.8l) . is a consequence of Theorem 4.6. of j5]. Uniqueness 
of a can found in [2]. Point (iii) gives a Primal-Dual optimality condition for fl3.7p 
and fl3.ip . it is a corollary of Proposition 13.31 it is very classic in the usual Monge 
case (see [1], Corollary 2.1.). The differentiability of u |cr|-a.e. is a well known 
result in the L^ theory of optimal transport, moreover it has been proved that (see 
for instance |25) for both properties): 

{z ejx, i/[, u{y) - u{x) = \y-x\) Du{z) = . 

Then by (i) and (iii), ^ = Du and so by Theorem 13.21 Du = D\„\u |cr|-a.e. □ 

Remark 3.5. In addition to point (iii), we have: 

(p, u) < j \y-x\ dj{x, y) + j g du^ 

for any admissible (m, 7 , z/) for (13. ip and (13.7p . Equality implies optimality of u 
and ( 7 , p). 

Let 7 be an optimal plan for (13.7p (and let u the corresponding boundary mea¬ 
sure) we can decompose 7 in four parts according to the origin and the destination 
of the mass transported 

7 = 7** + Ihi + 7*6 + Ibb 

with 7 jj = 7 |r 2 xo, 76 * = 7 |aoxo, lib = 7 |Ox 9 o, Ibb = 7 |aox 9 o, 
where i is for interior and b is for boundary. We have z/"*" = T^jiibi + Ibb) and 
p" = Tr^ilib + 766), then: 


\x-y\dj{x,y) -f / gdu^ = 


'QxQ 


Ian 


/ \x-y\d'jii + 

Iflxfl 

\x-y\d'yib+ / 


' 0x90 


/ {\x - y\ + g{x)}djbi 

'aoxo 

{\x-y\ + g{x)}d'ybb. 


'90x90 


Note that /gf 2 x 9 o'il^ “2/1+ 9 {.^)}dlbb{x, y) >0 and 7 — jbb is still admissible for 
(13.7p . Then since 7 is minimizing we have 


/ 90x90 


{\x-y\+ g{x)}d'jbb{x, y) = 0 . 
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From now on we will assume that ^hb = 0 (Otherwise spt( 76 ;,) C {(x,x) : g{x) = 
0}). The minimizing properties of 7 and the duality of Proposition 13.31 permits to 
characterize the points in the supports of 7 ;,* and 7 jf,. 

Proposition 3.6. We define the following multivalued maps: 

Pg (y) = {z ^ dfl : \y — z\+ g{z) < \y — x\ + g{x) \/x G 90}, 

p~{x) = {a; G dfl : \x — u\ < \x — y\ Wy E 90}. 

If (x, y) G spt 7 fei then x G Pg{y), if {x, y) G spt 7 jb then y G p~{x). 

Proof. We prove the statement about jhi, the other being similar. The part of the 
cost we look at is 

/ {\x-y\+g{x)} djbi{x,y). 

J dflxD. 

By dehnition of p'^ we have that \z — y\ + g{z) < \x — y\ + g{x) 7 f,j-a.e.(x, ?/) for 
all z G Pg{y). Then, for all measurable selection s of if we replace 'jbi by 
fbi = (s X id)]j^n‘^%i we obtain a new 7 admissible for fl3.7p and with lower cost. □ 

For later use we set: 

dg{y) := mm{^(x) + |a; - y\} = g{z) + \z - y\ Vz G Pg{y). 

The following Lemma will also be needed; 

Lemma 3.7. Let ( 7 , p) a couple of solutions of lid. 71) such that jbb = 0. If 
admits a solution such that u > 0 m O then: 

n is a non-negative measure and 7 = 7 ii + jbi 

with marginals p'^ = 71 ^ 7 , p~ = and n = 

Proof. Using Lemma [TH we get for 7 jb-a.e.(x, 1 /) the equality —u{x) = dQ{x,y). If 
ti > 0 this imply that 7 *;, = 0 and iy~ = 0. □ 

4. Some general results about the PDF 

The following result shows that, if it exists a solution to fll.ip . this solution is 
non-decreasing in time and unique. 

Proposition 4.1. Suppose that f^ > p and Uq > u^. Assume that for i = 1,2, 
{u^,Pi,Vi) are solutions of (PDF) satisfying U.Il) . (C), (Bl) with f = /*, initial 
and boundary conditions: 

u^{x) = 0 (z/*) -a.e.x, u^{x) = g{x) {n^)^-a.e.x, 

u\x, 0 ) = Uq Vx E Q. 

Then we have P > P. 
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It is a classical result when / is absolutely continuous, the following proof is 
adapted from jT2]. It requires the following lemma. 

Lemma 4.2. Let v : [0, T[xI7 — )■ M 6e such that v G L^{0, T, and 

dtv G L^{0,T, A4(Ll)). Then we have: 

^ dx^ = ‘^{dtv{x,t),v{x,t))M{n),cun) inV'(]0,T[). 

Proof. This is easily seen since applying Theorem 2.3.1 of |22|, it exists <I>„ G 
T>(fix]0,T[) such that: 

[ ||<h„(-,f) - n(-,t)||oodt-t 0, [ - dtv{-,t)\{n) dt ^ 0. 

Jo Jo 

□ 


Proof, (of Proposition 14.ip By Theorem 13.21 we have for a.e. t g]0,T[: 

uX-,t) G argmax{{fl - dtu\-,t),v)M(n),Ci{n) ■ v e Cipii^), 0 < v < g on 911}. 

Let u'^{x, t) = max{n^(a;, t), v?{x, t)} and u~{x, t) = min{n^(a;, t), u^{x, f)}. Using 
the optimality of leads (/^ — dtu^,u'^ — u^) < 0 a.e. t and as > p and 
we have: 

{P — dtP, u'^ — P) <0 a.e. t. (4.1) 

On the other side, we have: u~ — P = {P — n^)l{„i<n 2 } = {P — n’'')l{„i<„ 2 } = 
P — and a.e. t. These equalities imply 

{P — dtP, u~ — P) = {p — dtP^, P — up < 0 a.e. t. 

Combining this last inequality with fl4.ip and using Lemma [4.21 we get: 

- m^||^2 = {dtW^ - dtP.u-^-P) 

= {—p + dtu'^, «■*“ — u~) + {p — dtP,u~^ — P) < 0 a.e. t. 

So lln’*' — n ^||^2 is constant in time and as at time t = 0, it is zero, we get P > P 
at any time. □ 

We are now able to prove Theorem 11.11 As in j^H], we introduce the following 
functional on C6(0): 


T / \ _ j 0 if n G and 0 < n < on dfl, 

} -l-cxD otherwise. 

For any p G A1(0), denoting dloo the sub differential of Joo, by dehnition: 

p G diooiy) -v^ n G argmax{{p,w) : w G Cipi{fl) and 0 < w < g on 90}. 

For the record we write the following optimization problems which corresponds 
to fl3.ip . fl3.2p and fl3.7p with p = ft — dtu. We will make an intense use of 
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these new notations in the present and following sections. Recall that the equality 
niin fl4.3p = niax fl4.2p = niin fl4.4p holds with; 


— dtu{-,t),v) : v E 0 < v{x) < g{x) on dfl}, 


(4,2) 


mm 

crGM(0,K'i),i/GM(an) 


d|cr| + / g{x)di'^{x) : —diva = ft — dtu{-,t) — u in 

JdQ 


(4.3) 


mm 


7GAl+(OxO),i/GAl(9n) UqxQ 

We have the following result: 


x-y\dj{x,y)+ g{x)du^{x) : 
Jan 


vrjy = dtu + i/+, vrjy = /* + z/ } . (4.4) 


Theorem 4.3. Let u G T°°(]0, T[, hh^’°°(ff)) be such that dfU E h°°(]0, T[, A4(ff)) 
and 

\Du{x,t)\ < 1 a.e. {x,t), 0 < u{x,t) < g{x) on 9f2x]0,r[. 

1) (i) If {u, y, v) is a solution of U.l\} satisfying U.2\} then: 


ft - dtu{-,t) E dlcx,{u{-,t)) a.e. t e]0,T[. 


(ii) Assume u is non-negative, m(-,0) = 0, and ft — dtu{-,t) E dIoo{u{-,t)) 
a.e. t g]0,T[. Then it exists y E L°°(]0, r[, 2W’^(f2)) and z/ G 
L°°(]0, T[, Wl+(c)f2)) such that {u,y,u) is a solution of \1.1\) . 

2) Assume the conditions above are satisfied. 

• (Unigueness of u and y) The function u is unigue. Moreover, to each 
V corresponds a unigue y, and taking any jt such that ( 7 t,Pt) is a 
solution of iji4-4\ ) ^ following formula holds: 


{yt,T)M{n),cyn) ■= 



02 Jo 


(p((l - s)x + sy)\y - x\ dsd'yt{x,y) V(p G Cb{fl). 


For a.e. t g]0,T[, u{-,t) is space differentiable yt-almost everywhere 
and: Dfj,^u{x,t) = Du{x,t) yt-a.e. x. 


Remark 4.4. a) In the original article [29], in case g = d and / G T°°(]0, T[, L°°(f2))', 
L. Prigozhin, in the hrst place, proved (a similar result to) point 1). In the case 
he considered, y is expected to be in L°°(]0, T[, IT^’°°(ff))' and (PDE) has to be 
understood in the following sense; 

{dtU -f,ip) + {y, Du ■Dip)=0 V(p G L^QO, T[, L°°(f2)). (4.5) 

In our case, fl4.5p can be recovered by extending y in L°°(]0, T[, L°°(f2))' using the 
Hahn-Banach Theorem (cf Theorem 14.31 point 2)). 

b) The proof of l)(ii) below says actually more than required. Assume that ft — 
dtuf, t) E dIoo{u{-, t)) for a.e. t g] 0, T[ and take, for a.e. t, {at, nt) any solution of 
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(14.3p with Pt > 0. Then, under the assumption of the theorem, setting fit = \at\, 
we have that {u, fi, v) satishes (PDE) and 

t)| = 1 /i-a.e. (x, t), m(x, t) = g{x) p-a.e. (x, t). 

Moeover we have: at = 

Proof. l)(i) is an immediate consequence of Theorem 13.21 (since z/t > 0). 

1) (ii) With these assumptions, u{-,t) is optimal for (14.21) for a.e. t. Take {at,Pt) 
a solution of (14. 3 p with p non-negative: this is possible by lemma 13.41 and 13.71 By 
Theorem 13.21 fiii. setting fit = |cTt| we get the result. 

2) The uniqueness property of u comes from Proposition 14.11 while the properties 

of fi are contained in Lemma [TH as {at := D^^ufit,Pt) is a solution of (14. 3 p by 
Theorem 13.21 For the last property, see Lemma 13.41 □ 

5. The case of a finite number of sources 

In the spirit of |3] and [1], we are now looking at (dH]) when ft = 

As it is constant in time we will often write / for ft. 

In the next two lemmas we develop an heuristic of the shape of solutions for this 
special /. Starting from this idea we then show existence of a solution in Theorem 

El 

Assume for a while that a solution u G L°°{0,T,W^’°°{Q)) of (II.ip is known. 
By the previous section, we know u{-,t) is a solution of (14. 2 p for a.e. time, it’s 
non-negative and non-decreasing in time. Moreover take (yt, Pt) a solution of (14. 4 p 
for a.e. t G [0,r] with the following decomposition given by Lemma [3.71 

= 7ii,t + lbi,t and Pt = > 0. 

The Lemmas 15.11 and 15.21 will give us some clues to guess the shape of u: 

Lemma 5.1. Let us set: 

rj{-) ■= uivjr), u{x,t) := max{rj{t) - \x - yj\,0}. 

Thenu{-,t) is in Cipi{Ll) for all t g]0,T[, G C([0, T], M"*") and they satisfy: 

(i) u{x,t) > u{x,t) for all {x,t) with equality for x = yt and any x G 
spt(c)4M(-,f)) Uspt(pt), 

(ii) rj{t) G [0,d+(i/j)] with rj{t) = d+{yj) = g{x) -F |x - yf 'ybi,t-a.e. {x,yj) 
for a.e. t G [0,T]. 

Moreover 0 < ti(x, t) < g{x) on dVL for all t and u{-, t) is also optimal for ^.2^ . 

Proof. (i) The inequality comes from the Lipschitz property of u and the def¬ 
inition of rj. The equality at every yj follows from the dehnition rj{t) = 
u{yj,t) > 0. Then by Lemma [3.41 

u{yj,t) - u{x,t) = \x - yf 7 ra.e.(x, y^). (5.1) 

This gives the equality on spt{dtu{-,t)) U spt(r't) since = dtu{-, t) + Pt. 
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(ii) Consider x G p'^{yj), by the boundary condition (Bl) and the Lipschitz 
property of u: 

rj{t) := u{yj,t) < u{x,t) + |x - yj\ < g{x) + |a; - yj\ = for all j = 1, 

This implies u{x,t) < g{x) on dQ. Moreover, by Proposition 13.61 combined 
with fIS.ip every inequality above becomes an equality for i^^yj)- 

As is admissible for (14.21) and equals on spt(/) U spt{dtu{-,t)) U 

(spt(z/i)), it is optimal. □ 

For j = 1,..., k, let us introduce the following subset of hi: 

Aj{t):={xeQ: rj(t) - \x - yj\ = max{rn{t) - \x - yn\,0}}, (5.2) 

n 

k 

so that u{x,t) = - \x - yj\)lAj{t){x). (5.3) 

i=i 

Lemma 5.2. Assume for all j = 1, ...k, rj is derivable for almost every t g]0, T[. 
Then for all t > 0 and a.e. x: 

k 

dtu{x,t) = '^rj{t)lAj(t){x) a.e. {x,t)- 
i=i 

Proof. We have Int(Aj(t)) = {x G : rj{t) — \x — yj\ > Yna.-x.m^j{rm{t) — |x — 
ym|, 0}} and the boundary of Aj{t) is negligible. Indeed the sets Amif) D Aj{t) = 
{x G : \x — ym\ — \x — yf = rmif) — rj{t)} and {x G : |x — yf = rj{t)} are 
negligible as d > 2 and hi is convex. 

Let X in the interior of y4j(f), then rj(f) —|x—i/jl > maxm^j{rm(t) —|x —i/m|, 0}. As 
rj G C([0,T]), this inequality remains true at time t + h with \h\ small enough, as 
a consequence if is derivable on t: lim/i^o = liiu/i^o ^ 

rj{t). We then get dtu{x, t) = rj{t) for all x G Int(Aj(f)) for almost every t. 

□ 

In the sequel, we are going to see that, with the appropriate choice of rj, the 
application u = u is the unique solution of (II.lb . 

By Lemma l4.ll if u is a solution of (11.ip . the function rj := u{yj,-) are non¬ 
decreasing, moreover by the initial condition (I), rj{0) = 0. Take y* optimal for 
(14.4p . We deduce that, for small times rj{t) < d'^{yj) and by Lemma [5Tl for such 
t, 'ybi,t and Ut are 0. So that, ii u = u, then, for small t, the map has marginals 
Yl^j=irj{t)lAjit){x) a.e. (x,f) and / = X]j=i hy (1^ and (15^ . this implies: 

rjit)\Ajit)\ = Cj Wj = l,...,k, (5.4) 

k 

lt{x,y) = ^rj(f)lA,.(t)(x) 0 Sy.{y). 
i=i 
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Finally, every time (15.4j) holds, rj remains strictly non-decreasing and by Lemma 
15.11 when Vj reaches d'^{yj), it cannot increase anymore. 

All this remarks leads us to look at the ODE (15.5|) H. 

Lemma 5.3. There exist times {ti,...,tk) G and functions rj G C^{0,tj) 0 
C°([0,-|-cxo)) which satisfy: 

( o(^) = 1^/(4)! 

\ 0(0) = 0, (5.5) 

[ rj{t) = d+{yj) = min„ean{5'(a;) + \x - yf} Vt > tj 

with Aj{t) defined above in / (5.^l) . Moreover yj G Aj for all t and if t < tj then 

0<|A,(t)|. 


Proof. Step 1: Let mi = miuj^j \yi—yj\, m 2 = miuj d{yi, dQ) and m = min{mi, m 2 }, 
also consider c = max, c,. For small times i.e. for 


t < 


xJd 


,m. 


:= to 


{d+l)c 2 ' 

the functions ri{t) = (f^(d + l)t)^/'^'''^ are solutions of the ODE above with the 
correct initial data. 

Step 2; Starting from to the sets where ri{t) — \x — yi\ >0 may start to intersect 
or may touch the boundary of D and for the existence of solutions we appeal to a 
standard existence theorem. This requires that we prove that the functions 


are continuous B First we prove that, for every i G {1,..., fc}, (ri,... r^) i-)- 
|Aj(ri,... rk)\ is continuous on M"*" x • • • x M’*'. 

Let e = (ei,..., £fc). Dehne e = max^jej, e, — Sj} and e = miujle,, Si — Sj}. 


A(Li + ^1) • • • , Pfc + ^k) — 

= {x G D : 0 < Tj - |x - i/i|, Tj + Ej -\x- yf <ri + ei-\x- yi\ V j} 
C {x G D : 0 < Tj -f F — |x — ?/,|, rj — \x — yf < ri + e — \x — yi\ \/ j} 

= Aj(ri,... ,ri£,... ,rfc) 

and similarly Ai(ri,..., + e,..., r^) C Ai{ri + £ 1 ,..., + £fc). 

If \ 0 then Ai(ri,..., + e^+i, ■■■,rk) C Aj(ri,..., -1- .. ,rfc), and 

Ai{ri ,..., Tfc) = n„Ai(ri,..., ..., Xfc), it follows that 

\Ai{ri,...,rk)\ = \im\ Ai{ri,... ,ri + £n,... ,rk)\. 

n^O 


^ This ODE first appeared in [3], see also [4]. 
^This holds true if d > 1. 
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If Z' 0 then A(n,..., r* + , r^) C A(ri, • • •, + £„+i,..., r/,), and 

Z(li, ... ,rfc)\({x e : 0 = ri-\x-yi\}'OUj{x G : r —|a:-i/j| = ri-\x-yi\]) = 

Un Ai{ri, ... ,ri + ... ,rfc), 

and sinc^ 

|{a: G : 0 = r* — |a: — yi\}\ = 0 and |{x G : Vj — \x — yj\ = r* — |x — i/i|}| = 0 
it follows that 

\Ai{ri,...,rk)\ = \im\ Ai{ri,... ,ri + En,... ,rk)\. 

n—>-0 

Step 3 ; To get the existence of a solution of the ODE, it only remains to show 
that the measures of Ai{t) do not tend to 0 when rj(t) < d'^{yi). By the triangular 
inequality, if A* 7 ^ 0 then yi G Z- If |A(t)| = 0 then there exists at least 

one j ^ i such that — ri{t) = \yi — yj\. Indeed, if for all j 7 ^ z we 

have 0 > Lj = — ri{t) — \yi — yj\ then there exists 0 < I such that 

B{yi,l) C A('S) for s close enough to t which prevents the measure of Ai from 
going to 0 . 

Then we proceed by contradiction and assume it exist t and zq such that 
\imt^t\Aio{t)\ = 0 and ri^{t) < d^iyi^). Let I{t) := {i : ri{t) < d+(z/j)}. Since 
t HA ri{t) is non-decreasing 

5 ^ \Mt)\ > Y. ° 

SO there exist at least two indices {m,j) G /(^ x I{t) such that 

( 1 ) - rj{t) = \yj -ym\, 

( 2 ) limi^t \Aj{t) \ = 0, 

(3) limi^t \Am{t)\ = L > 0 . 

But this is impossible since, taking a derivative (as j G we also deduce 

— rj{t) = — 00 . Which concludes the proof. 

Step 4; Now we observe that 

D > 

|Z| " Z| 

and then reaches the value d'^{yi) in hnite time ti. 

□ 

We are now able to give the solutions of fll.ip and the associated optimization 
problems; 

• Set u{x,t) ■= ~ ~ yj\)^Aj(t){x) with Vj and Aj as in the 

previous lemma; 

^and here we are using both the convexity of LI and the dimension 2 < d 
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• Take Ii{t) = {j : t < tj} and l 2 {t) = {j : t > tj} two families of indices 
and set for any t ^ {tj : j = I,k}: 

Vt '■= Pj with Vj G M.'l{dVL) snch that: yj{dVt) = Cj, spt(pj) C PgiVj), 




I kin 

J6/i7) ' 

For a.e. f, define pt £ as: 


j&Ht) 


ip{x)dfit{x) = 



'02 Jo 

The following resnlt holds: 


</ 9 ((l - s)a; + sy) x \y - x\ dsdjtix, y), V(p G C 6 (hl). 


Theorem 5.4. With the above definitions, U.W holds, moreover: 

1) The triplet [u, p, v) is a solution of U.l\) . with uniqueness on u, 

2) For a.e. t G [0,T], u{-,t) is a solution of and ( 7 ^, vf) is a solution of 

U-4\ ), 

3) The couple ,f)ptWt) is a solution of a.e. t E [0,T]. 

Proof. Step 1: We show hrst that fll.2p is satished. Indeed for a.e. t g] 0,T[ by 
(ODE): 

0 < / dtu{-,t) = V / rj{t) dx = ^Cj< [ df{x), 

Jn j^j^ JAj j^j^ Jn 

by dehnition of u: 



and hnally: 


dpfix) = snp 


0=1 



02 Jo 


(p((l — s)x + sy) X \y — x\ dsdj^x, y) 


< diam(f 2 ) x 7 ( 0 "^) = diam(D) x / df{x) 


Jn 

Step 2: Note that, by Lemma [5.31 as yj G Aj{t), we have u{yj,t) = rj{t) for all 
time t g] 0,T[. Let ns now prove 2). Indeed, as for any j G /i, Cj = rj\Aj\, we 
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have for a.e. t g]0,T[: 


k k p 

~ *9t'u(-, t), m(-, f))_A4(o),ct,(o) ~ y ^ 1 “ y ^ / 

1 = 1 1 = 1 16/1 


u{x,t)rj{t) dx 


Y1 / ^ Cjrj(f) = / 1^ “ ^ (2^. 

ie/i -^^3 jeh jeh 16/2 


'02 


'02 


|x - 1/1 d'yii^t{x,y) + ^ 

16/2 ■ 


PgiVj) 


j&h j^h ^3 

\x — Hjl + g{x) diyj{x) 


x-y\ d'ya^tix, y) + \x - y\ d-fbi,t{^, y)+ 9{x) dvt{x) 


ivd 


I an 


\x - y\ djtix, y)+ g{x) dut{x). 


'an 


As we know u{-,t) is admissible by Lemma [5.II so. by duality (see Remark l3.5p . it 
is optimal and also is (yt, Ut). 

Step 3 : Lemma 13.41 gives a measure at G Adfe(r2,M'^) such that {at,iyt) is an 
optimal solution of fl4.3p for a.e. t g]0,T[: 


{(X, (f))M(nY,Ch{nY 



02 Jo 


0((1 - s)x + sy) ■ {y- x) dsd^t{,x,y). 


Then by Theorem 13.21 (7^ = D\^^\u{-,t)\at\ so 1) and 3) are proved (as yt = Wt\)- 


□ 


Note for later use the following estimations: 

Proposition 5.5. The following estimates hold for a.e. t g]0,T[.’ 

||M(-,t)||oo < llfi'lloo + diam(f2), ||T)n(-, f)||oo < 1, 

utid^l) < [ df{x), [ dtu{-,t)dx< [ df{x), 


dyt{x) < diam(r2) x / df{x). 


Moreover, we have: 


||5i//||L2(ox]o,r[) < ||M(•,^)||oo X / df{x). 

Jq 

Proof. The only remaining points are the estimates on ||n(-, t)||oo and ||c?iM||L2(ox]o,T[)- 
By the Lipschitz property of u and condition (Bl) taking x G hi and y G dO. leads: 

0 < m(x, f) < g{y) + \x — y\ for all t g]0, T[. This gives the first inequality. Let us 
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now show the last one: 

pT p pT k 

\\dtu\\L 2 inx]o,T[) = dtdx = / ^(rj(t))2|hl^(t)|rft 

J0 J O, j_^ J0 j—^ 

pT k k p 

= / = '^Cju{yj,T) < ||w(-,t)|loo / f{x)dx. 

Jo Jn 

□ 

6. A MORE GENERAL CASE 

Let /t = / G A^’''(f2), constant in time, we now aim to show that (II. Ih admits 
a solution (which will be unique by Proposition 14.Ij) . 

We approximate / by a sequence (/")n in such that; 

n 

r = nQ)<f{Q), r^finQ. 

i=l 

By the last section, it exists (m”, p"', in L°°(0, T; xL°°(0, T; Wl(f2)) x 

L°°{0,T] satisfying Theorem 15.41 and Proposition 15.51 

Proposition 6.1. (Convergence ofu'^, dtvP', u^) There exist u E L°°{0,T; VP^’°°(f2)) 
and p E L°°(0, T; Al(c?f2)) such that: 

i) u{-,t) ^ dfU E L°°{0,T-, JCl~^[dQ)) fl L^(f2x]0, r[) and u satisfies 

the conditions {Bi) and (/) of U.l\} : 

ii) up to a subsequence, m" converges to u for the strong topology o/L^(0, T, 
and dtu'^ converges to dtu for the weak star topology of L°°{0,T,JCI[Q)); 

iii) p{-,t) E a.e. t and converges to v for the weak star topology 

of L°°{0,T,A4{dfl)) (up to a subsequence) and {B 2 ) is satisfied. 

Proof. Step 1; By the estimates on {vT)m h is a bounded sequence in BV (f2x]0, T[), 
as a consequence (see Theorem 4 pi76 of |26)): 

—)■ M in L^(f2x]0,T[) and for a.e. t 

for some u E L^(flx]0,T[) and some subsequence {u'^'^)k of (ti"')n- Possibly ex¬ 
tracting again a subsequence, we can assume: 

u^'"{x,t) —)■ u{x,t) for almost every {x,t) E flx]0,T[. 

Step 2; By Ascoli-Arzela Theorem, for all t E [0, T[, («"'=(•, t))k admits a cluster 
point Vt E CipfiQ) (depending on t) for the uniform convergence. Using Step 1, we 
get u{x,t) = vfix) for almost every (x,t) G f2x]0,T[. This shows the uniqueness 
of the cluster point of ('u"'''(-, t))^ for the uniform convergence, so that: 

«"'=(•, t) —)■ u{-,f) Vt G [0,T[, uniformly in f2. 
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This gives (i?i) and (/). To get the convergence in L^(0, T,Cb{fi)) of we nse 

the dominated convergence Theorem and the following bonnd (n”'' is in Cip^{Q) 
and it is bonnded by || (71100 on dfl): 

snp t) — u{x, t)| < 2(diam(r2) + ||5'||oo)- 

t,X 

Step 3: L°°{0,T-,Cb{^)) being separable, (see for instance |22], Corollary 1.3.2. 
p 13), nsing the bonnd on dtu^, possibly passing to a snbseqnence, it admits a 
limit m G L°°(0, T, A^(r2)) for the weak star topology. It can easily be seen that 
m = dtu in the sense of distribntion. In the same way, dtu G L‘^{Qx]0,T[) and, 
np to a snbseqnence dtu^ dtU in L^. 


The point iii) is left to the reader. To get {B 2 ) jnst note that 

0 = lim g-Un) = {ly, g-u) 

n—>-+oo 

np to some snbseqnence. □ 

Remark 6.2. The proof shows also a different kind of convergence for {un)n'- up 
to a snbseqnence, converges nniformly to u for all t. As we will prove, u 

is the nniqne solntion of (II. ip . as a conseqnence the complete seqnence 
converges to u nniformly for all t. 


Proposition 6.3. (Convergence of and DgC) It exists p G L°°(0, T, Ad(r2)) 
and ^ G such that, up to a subsequence: 

g"' ^ g for the weak star topology 0 /L°°(0, T, Al(r2)), 

D^^u^g^ ^ f^g for the weak star topology of Ai{klx]0,T[Y, 


lim inf 

n—)-+oo 


dg^ > 


\^{x,t)\dg{x,t). 


'0 Jn 


'0 Jn 


Proof. The convergence of {g^n follows from its bonndedness. Moreover, we have: 


D^^u^{x,t)\^dg"'{x,t) = 


dg^<Tx df 


10 Jn 


'0 Jn 


so Lemma 3.3 pl3 of [6] applies and we get the rest of the proposition as: 


'0 Jn 


\Df,nu^(x,t)\dg'^{x,t) = 


dg^ 


>0 Jn 


□ 


Proposition 6.4. (Passing to the limit in the PDF) The following equality holds: 

g-a.e.{x,t), \Dfj^^u\ = 1 g-a.e.{x,t). 


Moreover: 


-div{Df,^ugt) = f - dtU - u 


m 


'x]o,r[. 
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Proof. Passing to the limit, it is easily seen that: 


-div(^/Zi) = f - dtu-o 


in 


^x]0,T[. 


So that is admissible for (14. 3 j) at almost every t g]0, T[. As, by Proposition 

16.11 is admissible for (I4.2j) at almost every t G]0,r[, by dnality we have: 

(/-atn(-,f),M(-,f));K(o)A(u) < supdi^D = infdiSD < / Ifi',t)\dnt+[ g{x)dvt{x) a. 


' an 


By Theorem 15.41 we have: 

[ {r - dt = 

Jo 

Then by Proposition 16.31 




'an 


( 6 , 1 ) 

g{x)di'"'{x, t). 


{f - dtu{-,t),u{-,t)) dt = lim / {P - dtu^{-,f),u^{-,t)) dt 

n^+oo Jr. 


= lim inf 

n^+cxD 


dp + 


'an 


g{x)du"'{x, t) > 


|(^| dp + 


'an 


g{x)di'{x, t) 


This implies that fl6.ll) is an equality and u is optimal for fl4.2p and v) is 
optimal for (14.3p . Then by Theorem 14.31 and Remark 14.41 we get the desired 
result. □ 

To conclude (recall Remark 16.2p , we have proved the following result: 

Theorem 6.5. Let f G A4’''(r2), the equations U.l\) admit a solution {u,g,u) 
satisfying Moreover: 

1) u is the unique solution of U.l\) . it is non decreasing. For a.e. t, the mea¬ 
sure dtu{-,t) is absolutely continuous with respect to the Lebesgue measure. 
Actually: dtu E L^(r2x]0, r[). 

2) z/ is non-negative, supported on p+(spt(/)). Moreover, to any u corresponds 
a unique g which can be built as: 

n 


{Pt,p) = 



Jo 


(p((l - s)x + sy)\y - x\ ds dpx,y) Mg E Cb(fl). 


where is any solution of ^.fD - 

3) For a.e. t g]0,T[, u{-,t) is space differentiable gt-almost everywhere and: 
Dfj,^u{x,t) = Du{x,t), gt-a.e. x. 

Moreover, for any time such that u{-, t) < g on dVL we have vt = 0, gt « 
Cn, and uniqueness on gt. 

4) Taking /"■ = Y17=ididyn any sequence converging (weak star) to f with 
yf'EQ and fn{kl) < f{Ll), u can be obtained as the uniform limit for all t 
of the sequence {u'^)n defined as: 

u'^{x,f) := max{r”(f) - dPx,y'f;),t)} 
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where rf G C^{0,ti) nC°([0,+oo)) satisfy the following ODE for some 
t'f > 0 ; 

f = ^ Vte]0,tf[, 

\ rr(o) = o; 

I = d+{y'^) Vt e [t^ +oo), 

A^{t) :={xen: - Ix-yfl = m^^{rf{f) - \x - y%t)}}. 

j 

Proof. It only remains to show the second part of 3). Indeed, when u{-,t) < 9, 
by Lemma [3.41 Ut = 0. So is the solution of the Monge transportation problem 
with marginals dtu{-,t) and /. As dtu{-,t) is absolutely continuous, so is pf (see 

IMlilSllillEQ]). □ 
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